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Abstract 



The role of the interaction in reaching and crossing the phantom 
divide line in the Chaplygin gas model is discussed. We obtain some 
necessary properties of the interaction that allow the model to arrive 
at or cross the phantom divide line. We show that these properties put 
some conditions on the ratio of dark matter to dark energy density in 
the present epoch. 

PACS: 98.80.-k, 98.80. Jk 

1 Introduction and preliminaries 

To describe the present accelerated expansion of the universe [T], many 
models have been introduced. In one of these pictures nearly 70% of the 
universe is filled by an exotic smooth energy component with negative pres- 
sure dubbed as dark energy. One of the candidates for the dark energy is 
the Chaplygin gas [2], a perfect fluid whose equation of state (EoS) is given 
by 

P d = -~. (1) 

Pd 

A is a positive real constant and and pd are the pressure and the energy 
density of the dark energy respectively in the rest frame of the fluid. The 
Chaplygin gas plays a dual role: it behaves as a dustlike matter in the 
early era (i.e. for small scale factor a), and as a cosmological constant 
at late times (i.e. for large values of a). However, this unified model of 
dark energy and dark matter, suffers from problems such as production 
of oscillations or exponential blowup of the matter power spectrum which 
is inconsistent with observation [3J. Besides this noninteracting model is 
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unable to describe the coincidence problem [3], i.e. this model cannot explain 
why the densities of the dark matter and dark energy are of the same order 
today. This lies on the fact that energy density of the non interacting cold 
dark matter (whose EoS parameter is w m ~ 0) redshifts faster than the 
energy density of the dark energy sector (whose EoS parameter satisfies 
Wd < —1/3), so the ratio of matter to the dark energy density is expected 
to decrease rapidly as the universe expands in noninteracting models. In 
addition this model is inconsistent with the phantom divide line crossing 
which based on astrophysical data seems to be occurred in the present era 

®- 

In this paper we assume that the universe is a spatially flat Friedmann- 
Robertson- Walker space-time which besides the dark energy component, 
whose EoS is given by ([1]), is composed of the cold dark matter interacting 
with the Chaplygin gas (which is taken as the dark energy sector) via the 
interaction term Q: 

p d + 3H( — 2 + l)p d = -Q 
Pd 

p m + 3Hp m = Q. (2) 

p m is the density of the pressureless (cold) dark matter. H is the Hubble 
parameter and "dot" denotes the derivative with respect to the comoving 
time. The presence of the interaction source, by permitting the energy ex- 
change between dark matter and dark energy, may alleviate the coincidence 
problem [6], meanwhile, as we will see, depending on the form of interaction, 
may allow the phantom divide line crossing. 
The Friedmann equations are 

—A 

H = -4tt( + Pd + Pm) 

Pd 

H 2 = ^{Pd + Pm). (3) 

o tt'2 

The condition p m + Pd = Pc '■= guarantees the flatness of the universe. 
The universe is accelerating if 

0<Pd< - pm + V ? m -±^. (4) 

In the absence of the interaction, i.e. Q = 0, p d is given by 

p d = VA + Ca-z, (5) 

where C is a real positive constant. The (EoS) parameter of dark energy is 

Pd A («\ 
w d := — = o • W 

Pd Pd 



2 



For ([5]) we have —l<Wd<0 for a 6 (0, oo), hence itf^ = —1 cannot be 
crossed in this model. As a consequence the EoS parameter of the universe 
w = QdWd, where Qd = ^ < 1, cannot cross the phantom divide line. Using 

lim a _ >00 pd = V^4, one can show that the ratio of dark matter to dark energy 
density, defined by 

Pm /„n 

r = — , (7) 

Pd 

becomes negligible at late time: lima^^r = 0. This was expected because, 
as was mentioned, cold dark matter redshifts faster than dark energy com- 
ponent. Now to investigate what happens in the presence of the interaction 
[7], let us consider the simple interaction term 

Q = 3TH Pd , T>0. (8) 

In this model the solution to the first equation in ([2]) is 



Pd = ViTr +c,a " 6(r+1) < (9) 

where C is a real positive constant. In this case in contrast to the noninter- 
acting case where Wd was restricted to the set (—1,0), in principle the EoS 
parameter of the Chaplygin gas can cross the phantom divide line (wd < — 1), 
for 

/ YA \6Trrry 
a > TTTTT^ ■ ( 10 ) 



C(l + T) 

Obtaining an exact solution for p m or H is not straightforward, but 
asymptotic solutions can be derived. Using 

dH 2 A 3H 2 . 

= _ 8yr + ) n 

ax pd 8-7T 



where x := In a, we arrive at lim^oo p c = v/QT + 1)A. In this limit the 
ratio of dark matter to dark energy, r := = ^ — 1, obeys lim a _ s>00 r = T, 
In addition the EoS parameter of the universe satisfies lim a _ s>00 w = — 1. In 
[8], where the role of the interaction on dynamical evolution of Chaplygin 
gas was studied in details, it was shown that this is a stable critical point 
corresponding to a stable attracting scaling solution in which although the 
Chaplygin gas crosses the phantom divide line but the universe will enter to 
a de Sitter phase and the big rip is avoided. 

So, the presence of the interaction term, which was first introduced to 
alleviate the coincidence problem, may allow the dark energy to cross the 
phantom divide line. In the next section we will find some necessary condi- 
tions for Q to allow the dark energy to reach or cross the phantom divide line 
and using an example show that there may be a deep relationship between 
the coincidence problem and Wd — — 1 in the present epoch. 
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2 The role of the interaction in crossing the phan- 
tom divide line and the coincidence problem 



In the presence of an arbitrary interaction, obtaining an exact solution to the 
Friedmann equations (to see whether the model can reach or cross Wd = — 1), 
may not be generally possible. In this part, instead of solving these equa- 
tions, we try to obtain some necessary conditions which must be satisfied if 
the model reach Wd = — 1. These conditions put some constraints on Q and 
the parameters of the model. Besides, depending of the form of Q, Wd may 
be a decreasing function at t = to, where to is defined by Wd(to) = — 1. In 
this situation w d crosses the phantom divide line. 

In the Chaplygin gas model, at t = to (where Wd(to) = —1), using ([2]) 
and © we obtain 

Pd(t ) = -Q(*o), p d (t ) = VA. (12) 

Note that nowadays (i.e. t ~ to) the dominate part of the universe is 
composed of dark sectors. In addition based on astrophysical data [I] : 
p d ^) = g ) therefore the second equation in ([2]) for a real Hubble parameter 
requires Pd(to) > 0. 

The time derivative of the EoS parameter, (JBJ), at t = to reduces to 

w d (to) = -2Ap-\to)Q{to) = -2Q(to)p d \to). (13) 

If Q(to) = then idd(to) = 0. In this case, from fl2J), we deduce Pd(to) = 
and to obtain the first nonzero time derivative of Wd at t = to we must 
continue this procedure: By taking another time derivative from ([6]), and 
then using p'd(to) = —Q(to), which follows from d5J when Wd = —1, Wd = —1 
and pd = 0, we obtain 

w d (t ) = -2p-\to)Q(to). (14) 

If again Q(to) = we get Wdito) = 0. For a differentiable Hubble pa- 
rameter, this method can be continued to give p^ = — Q^ n ~ l \ w d n \to) = 
2Ap d 3 (to) p d n \to) , where (n — 1) is the order of the first non zero derivative 
of Wd at to- Therefore finally we obtain 

w { d n \to) = -2p-\to)Q { - n - 1 \to). (15) 

If the dark energy sector reaches (or tends to) Wd = — 1 from the 
quintessence phase, depending on whether it will enter in phantom phase or 
not, one of the following predicates must be true (we assume that at to ; Wd 
is differentiable): 
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i) Wd tends asymptotically to —1. In this case, following (|15p . all the time 
derivatives of Wd vanish at Wd = — 1, hence Q and all of his time derivatives 
must also vanish in this limit: 

rl n O 

Q(*o) = 0, ~^{t ) =: QW(t ) = 0, Vn € TV. (16) 

For a derivable continuous this can only occur at to — > oo. E.g. the 
noninteracting Chaplygin gas model, Q = 0, belongs to this category and as 
we have seen Wd = — 1 occurs asymptotically. As another example consider 

Q = X(P d + Pd ) = H^d + l)p d , (17) 

where A is a constant. This is the interaction considered for some scalar 
field (cf>(t)) models of dark energy or inflaton, giving rise to the interaction 
X(j) 2 [9]. It is clear that at Wd = — 1, we have Q = 0. Using ([2]) one can show 
that the higher derivatives of Q also vanish at Wd = — 1, therefore Wd = — 1 
can only occur asymptotically. 

ii) The first nonzero derivative of Wd is negative and of odd order at 
Wd = —1, in this case after reaching, Wd crosses Wd = — 1 line and the 
Chaplygin gas enters the phantom phase. In this situation, Q{to) > or, if 
Q(to) = at to the first nonzero derivative of Q is positive and of even order 

Q(t ) > 0, or{Q(t ) = 0, and Q (2n) (t ) > 0, Q {k<2n) (t ) = 0}. (18) 

The example proposed in the previous section (jSJ), belongs to this cate- 
gory. The condition T > guarantees Q > 0. 

iii) Wd = at Wd = — 1, and the first non- vanishing derivative of Wd is 
positive and of even order . In this case the dark component does not enter 
in the phantom phase and Wd = — 1 is the global minimum of Wd- In this 
situation, with the same method used in ii), it is easy to show that at to 
we have Q = and the first nonzero derivative of Q is negative and of odd 
order 

Q(t ) = 0, Q( 2 " +1 )(t )<0, Q( fc<2 " +1 )(t ) = 0. (19) 

In this case the phantom phase is neither accessible for the dark energy 
component nor for the universe: w > —1. 

Therefore in the context of the Chaplygin gas model, reaching or crossing 
Wd = — 1 may not possible for an arbitrary Q, e. g. following (i) and (ii) and 
(iii) models with Q < do not arrive at Wd = — 1. Besides, via considering 
Q as a function of p m and pd, the above conditions give some relationship 
between the density of dark matter to density of dark energy at to- In this 
view one can determine or put some conditions on the value of r = — 
at the present epoch where based on some astrophysical data Wd ~ — 1 is 
happened. So there may be a deep relationship between coincidence problem 
and Wd = — 1 crossing observed in the present era. 
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To illustrate this point, let us consider the interaction |10j 

Q = H(\ mPm + \ dPd ), (20) 

and study the conditions required to reach Wd = — 1. Following (i) and (ii) 
and (iii) we must have Q > at to- F° r an expanding universe this implies 
A m r(t Q ) > -Ad- 
For 

r(to) = (21) 

Q = and Wd = 0, but using pd(to) = and p m (to) = — 3Hp m (to) we obtain 
Wd(to) = —3\dH 2 (to) and following (ii) (in order that reaching % = -1 be 
allowed) we must have A^ < 0. In this situation the dark energy component, 
remains in the quintessence phase. This possibility is excluded when one 
considers the interaction Q = XH(p m + p d ) which is adopted in some papers 

If \ m r > — Ad at to, then w d o cross the phantom divide line. In this case 
for A m > 0, and < we obtain a lower bound for r(to): 

r(t ) > (22) 

at transition time. 

Although obtaining an exact solution to the Friedmann equations in the 
presence of interaction ([20]) does not seem to be possible but let us examine 
whether these equations admit the series solution 

w d = -l + w d0 t a + O(t a+1 ), 
Pd = Pdo + Pdit 13 + 0(tP +1 ), 

Pm = PmO + PraXt 1 + 0(€ l+l ), 

H = Ho + Hrt + ^t 2 + 0(t 3 ), (23) 

at t = (we have taken to = 0). ([1]) leads to a = f3, pdo = VA and 

Pdi = W<K> 2 ; ■ We have also assumed that H(to) 7^ which lies on the 
assumption Cl d (to) 7^ 0. We first consider the solution with a > 2. Putting 
(|23p into the first equation in ([2]) and equating the terms with the same power 
of t, after some calculation we obtain X m Pm0 + ^dPdo = 0, and a — 1 = 7 
and subsequently by using the second equation in ([2|) we obtain a = 2 (viz. 
a > 2 is not allowed), 7 = 1, p m \ + 3H p m0 = and 2p d \ = -X m PmiH - 
By collecting these together, the series solution for a > 2 at t = (which is 
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also consistent with ([3])) is obtained as 

w d = -1 - 3X d H 2 t 2 + 0(t 3 
Pd 



s/A - ^X d H 2 VXt 2 + 0(t 3 ) 



P 



rn 



H = H + 4TT^VAt + O(t 2 ), (24) 

where the Hubble parameter at t = is determined to be Hq = ^ ^1 — \J~A. 
This solution corresponds to the system characterized by (|21|) . 

Note that p m \ + 3-ffoPmO = in this example leads to Q(io) = 0. This 
lies on the fact that a = 2 in (|23p. and therefore the Chaplygin gas remains 
in quintessence phase: w d (to) = 0. Hence (fT3j) results in Q(to) = w d (to) = 0. 
As we will show this is not true for other possible solutions. 

For a = 1 (corresponding to the aforementioned case A m r(to) > — A<f), 
one can similarly show that the following solution exists: 

w d = -l + ^=((X m -X d )y/A-^-X m )H t + O(t 2 ) 
V A V 8ir J 

Pd = VA+(^\ m -\ d )VA-^\ m ^H t + O(t 2 ) 

?1-Vtj + f(X m -X d + 3)VA - ^(A m + 3)) Hot + 0(t 2 ) 



H = H +(-^f+^VA\t + O(t 2 ). 



(25) 



According to the condition (ii), (I18p and subsequently ()25p are valid when 

X m H 2 > ^(A m - X d )VA. (26) 

This solution crosses uu d = — 1. If X m = we must have A^ > which is the 
same result discussed in the previous section. If A^ = 0, only models with 
A m > are allowed. If we adopt A m = A^ = A, then w d = — 1 crossing is 
possible whenever A > 0. 

After the dark energy component of the universe passes w d = —1, the 
EoS parameter of the universe may cross w = —1 too. Note that based on 
the values estimated for w d and Q d , it is also possible that the universe has 
an EoS parameter satisfying w = w d Q d ~ — 1 even in the present era, e.g. if 
we adopt £l d = 0.73 and w d ~ —1.33 [12], we obtain w = —0.97. 

For a general dark energy model it can be shown that 

= __LJii Q (27) 

w 3Hi-n d 3H Pc (i-n d y { ' 
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In the Chaplygin gas, from Wd^ld = w we derive 

(|27[) together with (|28p form a system (which is not analytically solvable) 
describing the behavior of w and £ld in terms of the scale factor. If the 
system arrives at w = — 1 at t = r from ()28p we obtain 

*M - If. ( 29 ) 

The universe enters in the phantom phase provided that w(r) < 0. By 
considering (j27j) . this implies 

Q > 3iJp c (l - fi d ), (30) 

at t = t. In terms of f2 m = ^ this inequality becomes 

Q > 3H Pm (31) 

at t = t. Putting ([20]) in this inequality gives 

(A m - 3)r > -X d , (32) 

which for X m > 3 and A^ < 0, gives a lower bound for r at transition time: 

r > - 3-^4, (33) 

while for A m < 3 and Xd > 0, — Decomes an upper bound for r. For 
X = Xd = X m , the upper bound becomes r < <p^. 

To study the behavior of the system near w = — 1 as before we may 
consider the following series solutions for H and O d at t = r (note that 
H(t) = 0), 

H = ho + h^t-rY + 0((t-r) u+1 ), u> 2, h x >0 
V d = u + Ul (t-Tr + O((t-rr +1 ). (34) 

/i and v are the orders of the first nonzero derivatives of H and £ld at 
t = t respectively. As we are studying the transition from quintessence to 
phantom phase we have taken h\ > 0. Using (|27|) and the equation 

2 H 

W = - 1 ~3H2> ( 35 ) 

we find that if /i 7^ 1 then we must have fj, = u (for details see |13j), but 
inserting the series solution in (|27p yields v = // + 1 which results in ^ = 2, 
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and [i = l. Note that v = 2 implies that the system cannot remain in 
quintessence phase. We also obtain ui = (3 — A m ) ho + (A TO — — 3) IiqUq. 
Using (|22]), ©, and © we arrive at 

# = ^o-^(^ I (3-A m )^ + (A m -A d -3)) (t-r f 
+0((t-rf) 

64tt 2 A ( tn . 64vr 2 „.\ . 

+0((i-r) 2 ). (36) 
/io is not determined, but h\ > gives 

647r 2 ^4 

(3-A m )^ < — — (3- A m + A d ). (37) 
y 

This inequality, which using r = 1 ~^ d can be shown to be the same as (|32p . 
determines the relation between the energy density at transition time and 
the parameters of the model. 



3 Conclusion 

We considered the interacting Chaplygin gas as the dark energy component 
of the universe. Some necessary conditions on the interaction, allowing the 
Chaplygin gas to reach the phantom divide line (wd = —1), were discussed. 
Using these condition we found that the ratio of dark matter to dark energy 
density, r, must satisfy some inequalities (or equalities) which can alleviate 
the coincidence problem, i.e. at = — 1, r is whether determined in terms 
of the parameters of the system or a nonzero lower bound can be found for it. 
Via some examples, we showed that how the form of interaction determines 
some of the behavior of the system such as crossing the phantom divide line 
or remaining in the quintessence phase. We also obtained a series solution 
to Friedmann equations at phantom divide line. At the end via an example 
we studied the possibility that the EoS of the universe crosses w = —1. 
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